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Distribution of localized states from fine analysis of electron spin resonance spectra of 
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We develop an analytical method for the processing of electron spin resonance (ESR) spectra. 
The goal is to obtain the distributions of trapped carriers over both their degree of locahzation 
and their binding energy in semiconductor crystals or films composed of regularly aligned organic 
molecules [Phys. Rev. Lett. 104, 056602 (2010)]. Our method has two steps. We first carry 
out a fine analysis of the shape of the ESR spectra due to the trapped carriers; this reveals the 
distribution of the trap density of the states over the degree of localization. This analysis is based 
on the reasonable assumption that the linewidth of the trapped carriers is predetermined by their 
degree of localization because of the hyperfine mechanism. We then transform the distribution over 
the degree of localization into a distribution over the binding energies. The transformation uses the 
relationships between the binding energies and the localization parameters of the trapped carriers. 
The particular relation for the system under study is obtained by the Holstein model for trapped 
polarons using a diagrammatic Monte Carlo analysis. We illustrate the application of the method 
to pentacene organic thin-film transistors. 

PACS numbers: 85.30.Tv, 73.20.At, 73.61. Ph, 76.30.-v 
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I. INTRODUCTION 

The electron spin resonance (ESR) technique offers a 
unique microscopic probe of carriers in semiconductors 
with unpaired spins. It measures transitions between the 
quantum levels rUg = ±1/2 in the presence of a magnetic 
field^"~. The spectrum of the transition between two 
quantum levels constitutes a J-function provided there 
is no external interference with the quantum levels of 
the spin system. However, the quantum states of carrier 
spins undergo a variety of interactions with the environ- 
ment. This interaction destroys the (5- functional form 
of the spectrum and broadens it. The broadening of the 
ESR signal is a result of two fundamentally different con- 
tributions from the medium. The first is a decay of the 
quantum levels caused by interaction with the excita- 
tions of the environment. This mechanism leads to the 
Lorcntzian shape of the spectral line. The second contri- 
bution is a result of the inhomogeneities of the medium. 
One example is the interaction with nuclear spins that is 
known as hyperfine interaction; here the inhomogeneities 
are caused by the probability distribution of the nuclear 
spin moments. In this case the spectroscopic signal is 
the sum of the contributions from spin systems located 
in different surroundings. The energy of signal in every 
surrounding is shifted by the local magnetic field depend- 
ing on the environment so that the summed signal has 
an inhomogeneous shape. 

As an example, electronic spin of cationic pentacene 
molecule isolated in a solution exhibits inhomogeneous 
broadening of the ESR signal that arises because of hy- 
perfine coupling with 14 proton nuclear spins. The sig- 
nal is constituted of a series of individual lines due to 



the hyperfine splitting. The envelope function of the sig- 
nal is roughly reproduced by a Gaussian^"^. This fea- 
ture is consequence of the central limit theorem (CLT). 
The local magnetic field for each electronic spin is caused 
by interaction with 14 proton nuclear spins and the ran- 
dom energy shifts of respective ESR signals are inevitably 
spread in accordance with Gaussian distributions, as a re- 
sult of the independent nature of respective nuclear spin 
orientations. 

We note that the anisotropic values of g-factors are av- 
eraged out (or narrowed motionally) by the rotation mo- 
tion of the molecules in solution. In contrast, solid-state 
organic molecular crystals exhibit ESR spectra composed 
of individual lines that are broadened due to the faster 
decay rate. Then, there might be the case that the indi- 
vidual lines have to become unresolved and the resulting 
ESR spectrum has to be very close to Gaussian envelope. 

Recent major developments in the field of organic thin- 
film transistors (TFTs) allow high-precision field-induced 
ESR measurement (which is referred in the following sim- 
ply as ESR) for the carriers in semiconductor crystals 
or films composed of regularly aligned organic molecules. 
In the measurements, carries are doped without introduc- 
ing any randomness by using the field-induced technique. 
The ESR signal in organic TFTs was first measured and 
analyzed in the groundbreaking study by Marumoto and 
his coauthor^^. The ESR signal observed in pentacene 
TFT at room temperature appeared to be narrower than 
that observed in solution. The authors claimed that 
the narrower linewidth should be an evidence of a spa- 
tial extension of wavefunction. According to CLT, the 
linewidth of a signal coming from charge distribution cov- 
ering N molecules is narrower by the factor l/yN. Then, 



assuming that the signal is Gaussian it was concluded 
that the wavcfunction is spread over A^ « 10 molecules. 

However, we note that the analyses were performed on 
the ESR spectrum with non-Gaussian lineshape. Subse- 
quent studies have shown that the field-induced ESR sig- 
nal and linewidth is temperature dependentiSr— . Typ- 
ical pentacene TFTs and rubrene single-crystal transis- 
tors exhibit sharp ESR signal whose single-Lorentzian 
linewidth presents motional narrowing^ effects with in- 
crease of the temperature. In case of the pentacene 
TFTs, the feature is well consistent with the ther- 
mally activated multiple trap-and-release (MTR) trans- 
port with the activation energy of about 10 meV in the 
high temperature range. In contrast, the narrowing ef- 
fect by the increase of temperature is not observed below 
around 50 K. Actually it was demonstrated by the con- 
tinuous wave saturation experiments that all carriers in 
the pentacene TFTs are localized at T < 50K and all re- 
laxation channels are frozen at such low temperatures^^. 
However the ESR spectrum is still deviated from the sim- 
ple Gaussian at low enough temperature. Therefore, the 
deviation of the ESR signal from Gaussian shape is not 
a result of relaxation or motional narrowing but should 
be associated with the nature of weakly-localized carrier 
states in the organic TFTs. It is of current important the- 
oretical and experimental challenges in materials physics 
to understand the carrier transport of organic electron- 
ics devices, as they permit productions of large-area and 
flexible electronic product o"^^'^^ . 

In the present paper we analyze the situation when 
the ESR signal of a semiconducting organic molecular 
system is a smooth curve that deviates from the Gaus- 
sian linewidth even at very low temperatures where the 
carriers are localized. Here we assume that peculiar line- 
shape is caused by a further specific inhomogeneity of the 
pentacene TFTs, as associated with the distribution of 
weakly-localized carrier states which are responsible for 
the device operation. Given this assumption, we have de- 
veloped a unique technique for obtaining the trap density 
of states from a few to tens of meVs in pentacene TFTsi^ 
the algorithm of which is described in detail in this pa- 
per. Note that the obtained energy resolution for the 
trap density of states is much higher than that by other 
methods based on transport or optical measurementii"— . 
The mathematical method suggested here is rather gen- 
eral and can be applied to analyses of broad spectrum in 
a variety of problems outside of those considered here. 

In Section [H] we study the deviation of the inhomo- 
geneous low-temperature ESR signal from the Gaussian 
shape. We show that the signal from the one pentacene 
molecule is very similar to a Gaussian showing almost no 
individual lines from hyperfine splittings. It appears that 
the ESR lineshapes of a signal from many independent 
traps of the same kind must be a Gaussian fsee lll AHIIBp 
whose width is uniquely determined by a single localiza- 
tion parameter N^ff, namely an effective number of sites 
where carrier is localized. Therefore, it is concluded that 
the non- Gaussian lineshape can be ascribed to the super- 



position of signals from different kinds of traps, where 
each kind of trap is described by its own localization pa- 
rameters Neff . In Section III CI we derive an explicit re- 
lationship between the shape of the experimental ESR 
signal and the distribution of the traps over the localiza- 
tion parameters N^ff. This relationship is a Fredholm 
integral equation of the first kind where unknown func- 
tion is a distribution of traps. Section IIIIl presents an al- 
gorithm to solve it based on the stochastic optimization 
method (SOM)24-27^ -y^e describe the SOM and present 
an analysis of its sensitivity to the experimental noise in 
Sections ImX] and [HTbI respectively. Sections HTClllD] 
present experimental details, handling realistic noisy ex- 
perimental data by SOM, and results for traps distribu- 
tions in pentacene TFTs. Section llH El presents methods 
that find the limits of the reliability of the distributions 
obtained. 



Section llVl shows how the distribution of the traps over 
the localization parameter iVg// can be mapped as a dis- 
tribution over the binding energies Eb- This transfor- 
mation can be formulated generally although a particu- 
lar implementation of the mapping requires the explicit 
Eb — ^ef f relationship between the binding energy Eb 
and the localization parameter N^ff. This relationship 
for a given model can be obtained using the exact numeric 
diagrammatic Monte Carlo method^^, the analytic mo- 
mentum average method^Si^, the coherent basis states 
method^i^, or numerous other methods (see [3j| for 
a review). We consider in Section IIV Al the model of 
two dimensional Holstein polaron in the field of an on- 
site attractive center. The distribution of the trapped 
states over the binding energies Eb in pentacene TFTs 
are shown in Section IIVBI Sections |V] presents a discus- 
sion of our results and Section IVII provides conclusive 
remarks. 



II. ESR SPECTRA OF TRAPPED CARRIERS 

IN ORGANIC SEMICONDUCTORS: 

FUNDAMENTAL KNOWLEDGE AND FURTHER 

GENERALIZATIONS 



In this section, we introduce the well-known charac- 
teristic features of the ESR spectra of a single molecule 
and a cluster containing several molecules (|II Al) . We then 
present an analysis of a noticeably different case where 
the carrier is localized on a single impurity in a crystal 
piBp . Finally, we consider the case of traps of different 
origin where we can introduce a relationship between the 
lineshape of an experimental ESR signal and the distri- 
bution of the impurities over different localization param- 
eters (|irc|) . 



ESR spectra for single molecule and a cluster 
containing several molecules 



In this section, we consider a molecular crystal in which 
single molecule contains so many nuclear spins that its 
ESR spectrum has an inhomogeneous Gaussian shape. 
A typical situation for a carrier trapped in a molecular 
crystal is that it is localized in a trap and its distribu- 
tion over the molecular crystal sites i is characterized by 
a probability distribution {pi}. The temperature is as- 
sumed to be sufficiently low that we can neglect the "ho- 
mogeneous" relaxation leading to the Lorentzian shape 
of the ESR signal. It is also low enough to avoid self- 
averaging of the inhomogeneities by the "motional nar- 
rowing" mechanism. 

In this case the lineshape of the ESR signal is deter- 
mined by the "inhomogeneous" broadening caused by the 
site dependent distribution of the hyperfine interactions. 
When the typical width of the individual spectral lines 
of the split with hyperfine interaction quantum levels is 
larger than the typical energy distance between these 
levels^>^, the lineshape of the ESR signal is Gaussian. 
This shape occurs when the carrier is localized in either a 
single molecule or a cluster containing several molecules. 

The case of a carrier trapped in a crystal is notice- 
ably different from the case of a cluster with several 
molecules. The probability distribution over N molecules 
i {pi, i = 1, N} is uniform pi — 1/N in a cluster. In con- 
trast, the probability distribution in a crystal trap pi, 
which is density of the carrier in given site i, is not uni- 
form, and the only restriction is the normalization con- 
dition J2iPi = 1- However, as shown below, the ESR 
signal of a carrier in a trap always retains the Gaussian 
shape and the width is uniquely determined by the car- 
rier probability distribution {pi}. 

The simplest ESR signal considered in our study is that 
for a single molecule. The fine structure of the ESR ab- 
sorption by a single molecule in a condensed environment 
is frequently blurred by the broadening of the hyperfine 
levels. The ESR signal from a single molecule in this case 
is Gaussian. The standard expression describing the hy- 
perfine structure of one molecule is^ 

ni/i riklk 

R{B) = J2 ■■■ Y. Pimi,.-.,mk)x 

mi=-ni/i mk=-nklk 

1 r 



''(^-EliA.m,) +r2 



(1) 



Here k is the number of groups of equivalent nuclei, rii 
is the number of equivalent nuclei in the ith group, li is 
the nuclear spin in the ith group, F is the linewidth of 
each peak, P is the intensity of each peak and B is the 
magnetic field. If protons (/ = 1/2) are the only param- 
agnetic nuclei, as is the case for pentacene molecules, P 
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FIG. 1: (color online) Simulated ESR spectrum (black solid 
lines) of 1, 2, and 1.54 pentacene molecules based on experi- 
mental hyperfine splitting of one pentacene molecule with fits 
by Gauss distribution (red dashed lines). 
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where C^'j "' ' are binomial coefficients. 

For the particular case of the pentacene molecule 
we set F — 0.02mT and use the coupling con- 
stants [Ai] i = 1, . . . , 4} and numbers of equivalent nu- 
clei {n^; i = 1, . . . ,4} reported in [Sj]. The ESR signal 
obtained from Eqs. ([1]) and ^ can be represented (see 
Fig. [T^) as a curve fluctuating around the Gaussian en- 
velope 



Go{B) = 
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2-Kal 
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(3) 



with standard deviation ctq = 0.554 mT. 

The standard situation, known from the physics of 
gases and solutions, is the case where the carrier is local- 
ized in a cluster containing N molecules and its density 
is spread over N molecules. In this case the signal retains 
its Gaussian shape with the width of the distribution re- 
duced by the factor N^^^. The hyperfine structure of the 
A'' molecules can be simulated by Eqs. ([1]) and ^ by 
replacing Ui — > Nrii and Ai — >■ Ai/N. Figure [Ud shows 
an example of the spectrum for N — 2 with standard 
deviation a ~ cro/%/2. It is clear that the oscillations 
around the Gaussian envelope are quickly suppressed as 
N increases. 

The shape and the 1/VN narrowing factor of the ESR 
signal for a carrier distributed over a cluster with N 
molecules follow from the CLT. The shifts of the sig- 
nal. Hi = {Bi — Bq) for the i-th molecule, is an indepen- 
dent random variable with Gaussian distribution R{y) = 
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FIG. 2: (color online) (a) Generated Gauss distribution with 
(T = 3, (b) distribution of a uniformly weighted pi = 1/9 sum 
of 9 normally distributed random variables, and (c) distribu- 
tion of exponentially weighted pi = exp(— 5i) {6 = 0.8) sum of 
9 normally distributed random variables. The squares show 
the distribution of the generated random variables and the 
solid lines are the fits by Gaussian functions. The envelope 
Gaussian function is normalized to unity in all panels. 



Go{y) with standard deviation (Tq. By the CLT, the dis- 
tribution R{y) of the random variable y = N^^ J2i=i Vi 
is also Gaussian with (Jn = (To/\/N. 



B. ESR spectra for trap in crystal 

The situation for a carrier localized in a trap in a crys- 
tal is different from the above situation with N molecules. 
The latter case assumes uniform charge distribution, and 
thus the CLT applies. In contrast, the distribution 
over molecules i in a trap {pi} is nonuniform with the 
probabilities pi subjected to the normalization condition 
^j^Pi = 1. Hence, we cannot assume a Gaussian line- 
shape for the ESR signal; the lineshape must be studied 
separately. 

Regardless of the lineshape, the probability distribu- 
tion {pi} unambiguously determines the linewidth of the 
ESR signal. The linewidth is characterized by the stan- 



dard deviation a which is the root square of the second 
moment of the linewidth. If we consider the standard de- 
viation of a signal in a trap a and compare it with that 
from a single molecule cto we can introduce the effective 
number of molecules Neff{{pi}) to describe the linewidth 
of the ESR signal from a carrier in a trap. The distri- 
bution of the ESR shift B is the same for each molecule 
i with mean (B) — Bq and variance CTq = ((-B — -Bo)^)- 
Since the hyperfine configuration of molecules are inde- 
pendent of each other the variables yi = {Bi — Bq) are 
independent for different molecules i. Hence, the stan- 
dard deviation cr({pi}) of the sum of random variables 
y = J2iPi{Bi — Bo) is related to the single-molecule stan- 
dard deviation ctq by the expression 



^(fe}) 
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/E^? 



(4) 



It is natural to define the effective number of molecules 
Neff{{pi}), corresponding to the charge distribution 
{pj, to be 



CT(fe}) - To/^N,ff{{p,}) 



(5) 



Then, the effective number of molecules Neff{{pi}) is 
unambiguously determined by the charge distribution in 
a trap {pi} 



Neff{{p^}) = 



Ep^ 



(6) 



To study the shape of the ESR signal for a trap with 
charge density {pi} we generated by a standard method^^ 
random variables {yi} following a Gaussian distribution 
with dispersion ctq = 3 (Fig. [2^). Then, we studied the 
distribution of the random variable y = '^^Piyi- By the 
CLT the uniform distribution pi = 1/iV leads to a Gaus- 
sian shape of the signal with the dispersion narrowed by 
the factor \/]V (Fig. [2}d). After performing simulations 
with a large set of different distributions pi we conclude 
that the distribution of the random variable y = ^^Piyi 
is always Gaussian (see Fig. [5J:) with -/Ve//({Pi}) defined 
by Eq. dH). 

Hence, we conclude that the shape of the ESR sig- 
nal for carriers localized in a set of identical independent 
traps is uniquely determined by the distribution density 
{pi}. It is always Gaussian with the standard deviation 
a defined by Eqs. (O and dH). 



C. ESR spectra for several kinds of traps 

Since the ESR signal for independent identical traps is 
always Gaussian, we assume that a non Gaussian shape 
originates from the superposition of the signals from dif- 
ferent kinds of traps. Indeed, each different kind of trap 
is characterized by a different probability distribution of 
the trapped carriers. 



To describe the experimental spectrum SexpiB) by the 
superposition of the ESR spectra G{B,^) for each trap 
type ^ we must choose a parameter ^ that unambigu- 
ously characterizes the spectrum G{B, (,). It follows from 
the analysis in Section Hi Al that the ESR spectrum from 
identical traps is Gaussian and can be characterized by 
a single parameter N^ff. This parameter reflects the 
spatial distribution {pi\ of a charge in a trap ([B]) and de- 
termines the narrowing of the Gaussian G{B,^ — N^fj) 
width a{{pi\) — dQ / ^J N ef f {{pi}) with respect to the 
ESR width CTo of a carrier localized on a single molecule. 
Therefore, the ESR signal for the same trap type, char- 
acterized by the same spatial extension ^ = N^ff , can be 
expressed as 



Here Bmin {Bmax) is the lower (upper) bound of a mag- 
netic field where the signal is larger than the experimen- 
tal noise. Xexp{B) are experimental data and X{B) is 
obtained from the distribution D{Neff) 



X{B) 



dG{B,N, 



ffi 



dB 



DiNeff)dNeff (11) 



G{B,Neff) 



'Neff 



{B~B„)^ 



(7) 



Introducing the distribution function D{N^jf) of the 
traps in pentacene over iVe// we can express the experi- 
mental signal £exp{B) in terms of the superposition 



^exp\^ ) 



G{B,N,ff)D{N,ff)dN,ff (8) 



which is considered to be a solution of the integral equa- 
tion. However, such naive approach leads to huge unre- 
alistic oscillations of the solution D{Neff). Instead, we 
need to apply one of the advanced techniques developed 
for such equations. 

Section HlI Al gives a general description of the stochas- 
tic optimization method (SOM)^iri22 for the solutoion of 
Eqs. ([S]) and ©. Section ITlI Bl applies the method to the 
analysis of the ESR data and demonstrates the influence 
of experimental noise on the reliability of the results. Sec- 
tion unci presents an experimental technique to obtain 
ESR spectra suitable for the fine analysis of the data. 
Section IIIIDI introduces an algorithm that implements 
the SOM and presents results for the trap distribution 
in pentacene TFTs. Finally, Section IIII El demonstrates 
the limits of the reliability of the distribution obtained 
by solving Eqs. dH) and ©. 



which is a convolution of the distribution function 
D{Nf,ff) of the traps and the Gaussian signal for a trap 
type characterized by the localization parameter A^e//- 

Most techniques for ESR measurements detect the 
derivative £exp{B) over the magnetic field B and hence 
the experimental signal Xexp{B) = d£ ^,j.p{B) / dB is re- 
lated to the distribution function of the traps D{Neff) 
via 



Xexp{B) 



dGjB.^Neff) 
dB 



Hence, to obtain the distribution D{Neff) we must 
solve one of the integral equations ^ and ©. The 
experimental signals Xexp{B) {£exp{B)) and the kernel 
dG{B,Neff)/dB {G{B,Neff)) are known functions and 
the distribution D{Neff) is to be determined. 



A. Method to solve inverse problem 

It is notoriously difficult to solve the Eqs. (jH]) and (|9]) 
because these equations belong to the class of " ill posed" 
problems. Naively, the true solution of the Eq. © 
D{Neff), being convoluted with the kernel, produces the 
function X{B), - which coincides with the given function 
B'{^eff)dNeff . (9) Xexp{B). Howcvcr, the general feature of the practical 

implementations of Eq. 1^ is that the knowledge about 
the function Xexp{B) is "noisy" and incomplete. Specif- 
ically, Xexp{B) is known for a particular set of points 
{Bi,i = 1,M} with some errorbars resulting from the 
experimental noise. In this case, to find a "solution", we 
can introduce the residual function 



A{t)^X{B,)-X,,p{B,), ^ = l,^f 



(12) 



III. FROM ESR SPECTRUM TO TRAP 

DISTRIBUTION: CONVERSION OF 

EXPERIMENTAL LINESHAPE INTO 

DISTRIBUTION OF TRAPS OVER DEGREE OF 

LOCALIZATION 

Equation ^ [^], where EexpiB) [Xexp{B)] and 
G{B,Neff) [dG{B,Neff)/dB] are known and D{Neff) 
is to be determined, is a Fredholm equation of the first 
kind. Naively, to find the solution D{Neff) we must 
maximize the inverse deviation 



and optimize a measure of the deviation of " the solution" 
X{B) from the given data Xexp{B). For example, we 
could maximize the inverse deviation (1101) 



M 



Ei^wi 



(13) 




^exp{B) - X{B) 



(10) 



Naturally, A(i) is never equal to zero at all points i — 
1,M when realistic noisy data {Xexp{Bi),i = 1,M} are 
considered. Hence, even the best measure Q™""^ is not 
equal to infinity. Therefore, the only feasible strategy for 
Eq. (jni) is to find a solution, that is "the best" in some 
sense. 



The above features are the characteristics of the class 
of " ih posed" problems for which we can not get an exact 
solution and can only find the "best" choice of D{Neff) 
for the given data set {Xexp{Bi),i = 1,M}. The naive 
approach, where we simply maximize measure (J13p . leads 
to unreasonable "solutions". Typically, they have huge 
fluctuations which exceed the true values of D{Neff) by 
several orders of magnitude. To get a reasonable de- 
scription of D{Neff) we must suppress this "saw tooth" 
instability. 

There are two different strategies. The first is the reg- 
ularization method, e.g. the popular maximal entropy 
method'^^ as one of many such methods^iS^. This ap- 
proach maximizes a measure that is similar to (|13p but 
modified in such a way that the solution is smooth enough 
to suppress " saw tooth noise" . The main drawback is 
that the solution is corrupted by the smoothing regu- 
larization procedure. The second strategy uses modern 
stochastic approaches to obtain many statistically inde- 
pendent solutions (see^ii^ and the references therein) 
whose linear combination smoothes the " saw tooth noise" 
without corrupting individual solutions. Since this ap- 
proach has been shown^I to be better for the solution of 
Eq. dni), we applied the SOM^i^^ as a particular example 
of the stochastic technique^. 

The function D{Nf.ff) is a distribution function and, 
hence, it is non negative D{Neff) > and normalized 



DiN,ff)dN,ff = I 



(14) 



to a certain number /. The noise present in experimental 
data {Xexp{Bi) , i — 1, M} makes the problem of normal- 
ization nontrivial. Some approaches to handle the nor- 
malization problem are suggested in the next chapters. 




FIG. 3: (color online) Test of SOM procedure to restore distri- 
bution of spatial extent of traps D{Neff) from "experimental" 
ESR data Xexp{B) ([9]) under different levels of noise. "Exper- 
imental" ESR data with noise are given in panels (al)-(dl) 
and actual (restored) spectrum D{x) is indicated by dashed 
(solid) line in corresponding panels (a)-(d). Noise level is 
(a,al) / = 0.0, (b,bl) / = 0.01 [s„ = 10"^], (c,cl) / = 0.03 
[s„ = 3 * 10"^], and (dl,d) / = 0.1 [s„ = 10"^]. Restored 
spectrum is obtained by solving Eq. (|9]) using SOM. Insets 
show details of corresponding spectra. 



B. Tests of SOM and role of noise in experimental 

data 

The SOM has been successfully applied to many inte- 
gral equations. The kernels of the equations are different 
from those in Eqs. ([5]) and (jH]). The exponential ker- 
nel K{y,x) = exp[— yx] was examined in [2j|37H55|, and 
various kernels ranging from the Fermi distribution to 
the Matsubara frequency representation are considered 
in (56i-59]. We test the applicability of the SOM to the 
kernels of Eqs. ^ and (|9]). We also show how the sta- 
tistical noise of the experimental data can obscure the 
information that can otherwise be obtained by solving 
these equations. 

To verify that the SOM is applicable to Eqs (|8]) and 
© with the kernel defined by Q we introduced a nor- 
malized to unity function D{Neff) (see the dotted line 
in Figs. 13] and H]) and generated a set of "experimen- 
tal" data {XexpiBi),i = 1,200} using relations (09]). 

Then, we attempted to find D{Neff) by solving Eqs. 
(|S]) and (jH]). For ideal data without noise in the set 



{Xexp{Bi),i — 1,200} we were able to restore D{Neff) 
successfully (see Fig. [3lja) and IUJa)). We did not find 
a significant difference between the results obtained by 
solving Eqs. dH) and dH). 

In spite of the extensive usage of SOM^^*^"— there 
is no precise understanding of how the errorbars of the 
data {Xexp{Bi),i = 1, M} obscure the solution D{Neff)- 
The difficulty is the dependence of the result D{Neff) on 
the level of statistical noise, the shape of D{Neff), the 
kernel, and even the number of points M in the data set 
{Xexp{Bi),i — l,M}. Therefore, below we present some 
examples showing how the noise of the experimental data 
obscures the solution of Eq. ^. These ad hoc examples 
indicate the general trends but must not be treated as 
a quantitative analysis of the influence of noise on the 
reliability of the SOM results. These data should not be 
used for the analysis of other cases. In Section IIIIEI we 
present a procedure to check the reliability of particular 
solutions. 

To introduce noise we used a sequence of random num- 
bers TZ(i) uniformly distributed in the range [-1,1] and 




FIG. 4: (color online) See figure caption in Fig. [S] 



generated data sets {XexpiBi) + {f/2)n{i),i = 1,200}, 
referred in the following as "experimental" ESR data, 
with different amplitudes of the noise /. The signal- 
to-noise ratio is defined as a ratio of the amplitude / 
and maximal absolute value MAX{| X^xpiBi) \} of the 
given signal Xexp{Bi): as s„ = f/MAX{\ Xexp{,Bi) |. 
The results presented in Figs [3] and S] illustrate the gen- 
eral trends. An increase in the signal-to-noise ratio s„ 
corrupts the solution for large values of A^e// first: the 
shape of the high-energy peak is not reproduced but its 
position is still correct. At higher values of s„ the shape 
of the low-energy peak is not reproduced. A comparison 
of Figs 13] and S] shows that the spectrum with a sharp 
feature at small A^e// (Fig. |3]) is more robust to experi- 
mental noise than that with a broad feature at small A^e// 
(Fig. 2]). Note that although the shapes of the high- and 
low-energy peaks are not reproduced, their positions are 
still approximately correct even for large values of the 
signal-to-noise ratio s„. 



C. Experimental data for analysis 

To conduct reliable spectral analysis as discussed 
above, we need high-precision ESR spectra for the car- 
riers in organic TFTs. We acquired the spectra by the 
following procedures: We used a commercially- available 
X-band (9 GHz) ESR apparatus (JES-FA200, JEOL) 
equipped with a high Q cylindrical cavity (Q factor 4000- 



6000 for the TEqh mode). We fabricated bottom-gate, 
top-contact pentacene TFTs with high mobility that are 
suitable for high-precision measurements. The device is 
composed of a 100-/im-thick poly(ethylene naphthalate) 
(PEN) film as a nonmagnetic substrate, a l-/im-thick 
Parylene C film as a gate dielectric layer (4.5 nF/cm^), 
and a 50-nm-thick pentacene film as the semiconducting 
layer. The gate, source, and drain electrodes are com- 
posed of vacuum-deposited gold films with a thickness of 
30 nm; this is much thinner than the skin depth of gold 
(about 790 nm). 

Since the field-induced carrier is accumulated only at 
the semiconductor/insulator interface, the ESR signal is 
proportional to the total channel area of the TFTs. We 
used a device with a width of 2.5 mm and a length of 
20 mm, the dimension of which is limited by the inner 
diameter of the ESR tube and the cavity size. We used 
a stack of ten sheets of TFTs for the high-precision ESR 
measurement, to obtain field- induced carriers ten times 
as large as those in a single sheet. The total carrier num- 
ber at Vg = -200V is estimated as 2.8x10^3. 

The semiconducting pentacene layer is composed of 
a uniaxially-oriented polycrystalline film where all the 
component pentacene molecules aligned with the molecu- 
lar long axes are roughly perpendicular to the film plane. 
In the measurement, a static magnetic field was set per- 
pendicular to the film plane to eliminate the anisotropic 
effect of the g tensor. A continuous-fiow cryostat was 
used for the low-temperature measurements. In the 
FESR measurements at low temperature, we first applied 
the gate voltages at room temperature (with the source 
and drain shorted) and then cooled the device to the set 
temperature, to avoid a delay in the charge accumulation. 
The temperature was stabilized carefully so that the fluc- 
tuation at 20 K was about 0.01 K, which minimized the 
effect of temperature-dependent spin susceptibility. 



D. Practical implementation of method: 
Distribution of traps in pentacene TFT 

Equation (jH]) is preferable for the practical impli- 
cation of the SOM. The problem with realistic noisy 
data from ESR experiments is that there is some un- 
certainty in the normalization and background of the 
experimental data. In the ideal case, implying the 
normalization of the distribution D{Neff) to unity 
and the conditions /_ G{B,Neff)dNeff — 1 and 

J_ dNi^ff J_ dz[dG{z^Nf.ff)/ dz\ = 1, we must nor- 
malize the experimental data as 



dB£,xp{B) = 1 
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FIG. 5: (color online) (a) Dependence of the fit quality Q™"^ 
and position of low-energy peak S(N — Nsff) on normalization 
/ of the spectral function. 



respectively. However, the normalization of the experi- 
mental data is not exact because experimental noise can 
lead to uncertainty of a few percentage points. On the 
other hand, the solution of Eq. ([9]) is sensitive to possible 
normalization errors. 

There is also a problem with the background. If the 
experimental data are obtained in the form £f,xp{B) (see 
e.g. Fig. [2]), there is no unique procedure to determine the 
constant background level that must be subtracted from 
the data to get a pure signal for the ESR transition. This 
uncertainty also increases the uncertainty relating to the 
normalization of the data. However, the problem of the 
unknown background disappears when Eq. ([9]) is used 
as the integral equation because lims-s-itoo '^exp{B) — 
for any constant background. Therefore, to analyze the 
ESR data for pentacene TFTs we used ^ where the only 
uncertainty is that of normalization. 

To handle the normalization uncertainty we can change 
the normalization in either Eq. (fT4)) or Eq. (|16l) . The 
two choices are equivalent because of the linearity of Eq. 
(P|). In practice we normalized the experimental signal 
as shown in (IT51) and varied the normalization / of the 
distribution density D{Neff) ([T^. We handled the nor- 
malization uncertainty for the result shown in Fig. [6] as 
follows. The spectrum at gate voltage -200 V was con- 
sidered (the result for this gate voltage is shown in Fig. [5] 
by solid line). The integral equation (jH]) was solved for 
different normalizations / and the solution with the nor- 
malization having the best deviation measure Q™"^ was 
chosen. 

Figure [S] shows the best inverse deviation Q™"^ (|10ll3p 
and the position of the sharp peak in D{Neff) versus the 
normalization /. It can be seen that the position of the 
sharp peak is sensitive to the value of the normalization 
/ and this may be a source of the volatility in the fine 



FIG. 6: (color online) Distribution of trap states in pentacene 
TFT versus spatial extent Neff of charge distribution in traps 
obtained from ESR spectrum of pentacene TFT at 20 K with 
gate voltage -200 V (solid line), -120 V (dashed hne), and -40 
V (dotted line). 



analysis of the ESR data. However, it can be shown 
that the suggested approach to the the normalization /, 
which leads to the best deviation measure Q™"^, is ro- 
bust and produces stable results. We have demonstrated 
this via analysis of the ESR data at different gate volt- 
ages (Fig. [nj. We found that the best normalizations are 
different at different gate voltages: I{V — —200) = 1.01, 
I{V = -120) = 1.038, and I{V = -40) = 1.03, re- 
spectively. However, the position of the sharp peak at 
low values of N^/f does not depend on the voltage V if 
at each voltage V we use the normalization I{V) that 
corresponds to the best deviation measures Q™"^ . Phys- 
ically, the sharp peak at low values of Neff corresponds 
to deep impurity levels that depend only slightly on the 
gate voltage. Therefore, its independence on the gate 
voltage in the fine analysis of the ESR spectra indicates 
the high stability of the procedure based on the suggested 
approach. 



E. Reliability of trap distribution result 

Since solving Eq. ^ is an "ill-posed" problem it is 
useful to understand how much information we can get 
from the analysis and to check how many details of the 
resulting distribution of impurities are reliable. The reli- 
ability can be analyzed by plotting the residual function 

The spectrum D{Nf,ff) (Fig. [5]) obtained by solving 
Eq. ^ has three peaks. Figure. [7^ shows the fit of the 
ESR signal using the distribution D{Neff) in Fig. [S) It 
also shows the separate contributions of the A, B, and C 
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FIG. 7: (color online) (a) Experimental signal (squares), 
fit by spectrum in Fig. [6] at the gate voltage -200 V (solid 
lines) and contributions from (see Fig. [6]) A, B, and C com- 
ponents (dashed lines). Residuals (|12p Xexp{B) — X{B) 
from (b) unbiased distribution (Fig. (Gjl, (c) best fit by one 
Lorentzian, (d) best fit by two 5-functions 0.68S{N — 4.5) 
and 0.325{N — 20), and (e) best fit by three (5-functions 
0.31S{N - 1.4), 0.515{N - 7.5), and 0.18S{N - 25.0). 



components of the distribution. To clarify which features 
of the distribution D{Neff) are reliable for the given level 
of noise in the experimental data we studied the residuals 
(IT2t Xe^cpiB) - X{B) (Figs.Eb-e). We can see that the 
quality of the fit by the SOM (Fig. Eb) is much better 
than that obtained by e.g., the Lorentzian (Fig. [T^) and 
two (5-functions (Fig. [TJi). The fit from three (5-functions 
gives a residual function (Fig. [TJb) as good as that ob- 
tained from D{Neff) in Fig. [51 Therefore, we conclude 
that, within the limits of the noise of the experimental 
data, the existence of at least three kinds of traps is a 
reliable result. 

We note that the distribution over the parameter Neff 
in Fig. [S] is free from any assumption about the shape 
of the distribution. Indeed, because of the noise of the 
current experimental data, the only reliable conclusion is 
the statement about the existence of at least three types 
of traps. However, a data analysis with less noise could. 



in principle, reveal more fine structure in the distribution 
function D{N^ff). 



IV. TRANSFORMATION FROM SPATIAL 
DISTRIBUTION TO ENERGY DISTRIBUTION 

In this section we discuss finding the distribution of 
the traps Z{Eb) over the binding energies Eb given 
the distribution D{Neff) over the localization parame- 
ter Neff. This transformation is trivial when there is a 
priori knowledge of the functional dependence 



Neff = N,ff{EB) 

Indeed, there is a balance relation 
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The above generic relation (|17l) can be obtained as a 
parametric function provided we can to choose an ap- 
propriate model and calculate N^ff = Neff{{V}) and 
Eb = EbUV}) in some domain {V} of the parameters 
V. 

The general strategy for the transformation D ^ Z 
does not depend on the model. However, the specific 
form of relation (fT7|) depends on the model chosen to 
describe the localization of a carrier in a trap. To study 
carrier trapping in pentacene TFTs we consider a model 
of a two-dimensional (2D) Holstein polaron in a field of 
an attractive center. This model is considered in Section 
IIV Ai and the results of the transformation for pentacene 
TFTs are given in Section HvB 



A. Model for traps: 2D Holstein polaron in field of 
on-site attractive center 

It is already well known that the behavior of a car- 
rier in pentacene TFTs can be described by a parti- 
cle in a system with attractive impuritiesiS-. It is also 
known that the carrier is subject to the electron-phonon 
interaction^. To model the behavior of a carrier in pen- 
tacene TFTs we chose the simplest Hamiltonian describ- 
ing a 2D Holstein polaron in a field of on-site attractive 
center 



H = -t J2(u) qcj + wph Ei Hh 



-lJ:^{bl+bi)4ci-Ucico 



(19) 
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FIG. 8: (color online) Dependence of (a) binding energy Eb 
in units of bandwidth Wb = 8f and (b) effective number of 
molecules Neff on absolute value of potential of attractive 
impurity | [/ | in units of Wb- (c) Dependence of effective 
localization number N^ff on binding energy Eb in units of 
bandwidth Wb- Curves are presented for A = (squares), 
A = 0.15 (circles), A = 0.5 (triangles pointing up), and A = 1 
(triangles pointing down). Lines are to guide the eye. Inset 
in (c) shows schematically model represented by Hamiltonian 



Here, c] {b]) is the creation operator for the carrier 
(phonon) in the i-th molecule. U is the attractive im- 
purity potential for the carrier cj at site and Wph is the 
frequency of the dispersionless phonon. The amplitude t 
describes the electron transfer ex t between nearest neigh- 
bor sites and the local Holstein coupling to the phonons 
is oc 7. The dimensionless electron-phonon coupling con- 
stant A is defined to be A = 7^/(4ia;ph). It is clear that 
for the chosen model the parameter set to determine re- 
lation (IT7|) is T' = {U,X} including the potential of the 
attracting trap U and the strength of the electron-phonon 
coupling A. 

To calculate the values of Eb = Eb{U, A) and N^ff = 
Neff{U, A) we used the direct space diagrammatic Monte 
Carlo (DSDMC) technique^^. Similar data can be can 
be obtained by the inhomogeneous momentum aver- 
age approximation method— 1^ and the coherent basis 
states method^^'^^. The data for Eb = Eb{U,\) and 
Neff = Neff{U, A) are presented in Fig. |5^ and[Ht>. The 
values of N^ff were determined by relation ^ from the 
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FIG. 9: (color online) Charge distributions around attractive 
impurity with potential U at A = 0.15 along [100] and [111] 
directions: U/Wb = -0.206 [EB/Wt = 0.0032, N^ff = 44.3] 
(circles); U/Wb = -0.25 [Eb/Wi, = 0.0113, iVe// = 11.2] 
(rhombus); U/Wb = -0.281 [EB/Wt = 0.0213, N^ff = 6.05] 
(squares); U/Wb = -0.375 [EB/Wt = 0.0703, iVe// = 2.38] 
(triangles pointing up); U/Wb = -0.59 [Eb/W^ = 0.2397, 
N^ff — 1.35] (triangles pointing down). 



charge distribution in the trap which was calculated by 
the DSDMC technique (see Fig.E]). 

To determine unambiguously the functional depen- 
dence Neff — Neff{EB) from the calculated relations 
Eb = Eb{U,X) and Neff = iVe//(C/,A) (Figs. [8ji and 
[Hb) we must decide which of the two parameters, A and 
U, is fixed and which is responsible for the variation of 
the physical parameters of the traps: the binding en- 
ergy Eb and the localization parameter Neff- A proper 
choice of the parameter responsible for the variation in 
the physical properties of the traps determines relation 
(|17p fully and unambiguously. Since the thin film in pen- 
tacene TFTs consists solely of pentacene molecules it is 
natural to assume that the value of A is one and the same 
for the entire film and the spread of the physical parame- 
ters of the traps is caused by trapping potentials of differ- 
ent origins. The trapping potentials of different origins, 
in turn, can be characterized by different strengths of the 
attractive potentials U- 



B. Energy distribution of traps in pentacene TFTs 

To analyze the ESR spectrum of pentacene TFT we 
used the electron-phonon coupling constant A = 0.15, es- 
timated from optical absorption experiments^*', and the 
hopping amplitude t = O.leV, obtained from band struc- 
ture calculations^ii^. Figure [TU] shows the distribution 
Z[Eb) of the trapped carriers over the binding energies 
in TFT at the gate voltage -200 V. We fixed A = 0.15, 
considered the dependence of Eb (bold line in Fig. [8^) 
and Neff (bold line in Fig.[5]D) on the attractive potential 
U, and obtained Neff[EB) (bold line in Fig. [St). Then, 
we used transformation (|18p . 

We found that the two discrete trap levels (A and B) 
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FIG. 10: (color online) Distribution of trap states in pen- 
tacene TFTs as function of trap energy Eb for gate voltage 
-200 V. 



peak at 140 ± 50 meV and 22 ± 3 meV, respectively. 
The broad feature (C) at the gate voltage -200 V is 
distributed between 5 and 15 meV, as presented in 
Fig. 1101 The low-energy profile prompts the anticipa- 
tion of tail states extending from just below the band 
edge, as has been discussed for amorphous semiconduc- 
tors, while the states are partially occupied up to the 
Fermi level at around Eb = 5 meV. These results are 
roughly consistent with the small activation energy of 
about 15 meV for the motional narrowing observed in 
[14|. and also with the potential fiuctuations by atomic- 
force-microscope potentiometry^'^. We note that the dis- 
tribution Z{Eb) gives relatively correct position of the 
trap levels, although the absolute value of the binding 
energy is rather model-dependent. 



V. DISCUSSION 

Weakly-localized in-gap states are expected to play 
crucial roles in the intrinsic charge transport along semi- 
conductor/gate dielectric interfaces in organic transis- 
tors. In practice, temperature-independent mobility is 
often observed in devices with high mobility and highly- 
ordered molecular interfaces, which indicates that the 
Fermi energy is just below the band edge^'*. 

To date, various experimental techniques have been 
used to investigate the interfacial trap density, such as 
deep-level transient spectroscopy (DLTS)^^, photocur- 
rent yield^, gate-bias stress^, and thermally-stimulated 
current^ experiments. However, the measurements are 
based on the charge transport, and it is strongly affected 
by the " extrinsic" potential barriers at grain boundaries 
and/or channel/electrode interfaces. In striking contrast, 
our method has a crucial advantage in its ability to dis- 



FIG. 11: (color online) Comparison of distributions of trap 
states in pentacene TFTs as function of trap energy Eb for 
gate voltage -200 V (solid line), -120 V (dashed line), and -40 
V (dotted line). 



close the spatial and energy distribution of shallow traps 
down to a few meV, based on a unique microscopic probe 
using electronic spins. In addition, the g tensor can be 
used to identify the molecular species around the trap 
sites. For the three types (A, B, and C) of trap states 
obtained, the g tensor should be common, considering 
the highly symmetric nature of the ESR spectra. This 
indicates that the trap states are extended over inherent 
pentacene molecules of regular orientation^iiS.. Of these, 
the deep discrete trap level (A) might be attributed to 
structural defects such as grain boundaries^. 

The shallow discrete level (B) and the broad fea- 
ture (C) might be ascribed to small defects such as 
molecular sliding along the long axes of the molecules^^, 
disorder induced by random dipoles in the amor- 
phous gate dielectrics2ii2^, thermal off-diagonal elec- 
tronic disorderiSiZij and the fiuctuations of the band 
edge2^. Notice that the regular orientations of the 
molecules are retained in the trap states as stated above. 
Although these assignments are rather speculative, we 
believe that further microscopic investigations based on 
this study will provide a comprehensive view of the 
weakly-localized in-gap states in organic transistors. 

The lack of sensitivity of the strongly localized states 
with small N^ff to the gate voltage, which is obvious 
from the physical point of view, is an indication of the sta- 
bility and reliability of the fine analysis of the ESR spec- 
tra. For shallow states with large iVe// and small Eb , an 
increase in the gate voltage adds low-energy states that 
participate in creation of ESR signal (Fig. [TT|) . This shift 
of the border where states are visible to the ESR probe 
indicates that these states are filled as the gate voltage 
increases. Hence, this behavior can be interpreted as a 
movement of the Fermi level which tends to zero energy 
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as the bias voltage increases. 

It is important to mention that the sharp peak of 
D{Neff) at Neff = 1.54 does not contradict the assump- 
tions used to derive the integral equations (|H]) and (0). 
The very essence of these equations implies that the con- 
tribution from each state with a given Neff is a Gaussian 
ESR signal. On the other hand, the signal at small N^ff 
is a more complicated function with fine features (see 
e.g. Fig. [1^ for N^ff = 1). The results for distribution 
D{Nf.ff) at small N^ff can be unreliable. However, as 
can be seen in Fig. [TJ:, the ESR signal for reasonable 
parameters is close to Gaussian even at Neff = 1.54. 
Therefore, the results for the distribution of the traps 
D{Neff) are valid even for small values of the localiza- 
tion parameter Neff . 



VI. CONCLUSIONS 

We have presented an unbiased analytical method for 
the processing of high-precision electron spin resonance 
(ESR) spectra, which allows us to obtain the distribu- 
tion of trapped carriers over the degree of localization 
and the binding energy. The first step is a fine analy- 
sis of the shape of the ESR spectra by the SOM, which 
allows us to split the spectrum into multiple Gaussian 
components each of which corresponds to a different spa- 
tial extension of the trapped carriers. The second step 
is the transformation of the distribution over the degree 
of localization into a distribution over the binding ener- 
gies via a system-dependent relation between the binding 
energies and the localization parameters of the trapped 
carriers. We have presented and discussed the fundamen- 
tal bases of the spectral analysis, detailed algorithms for 
practical applications, and discussed the robustness of 
the analysis to experimental noise. Although the method 
can be applied to many systems, we consider that it is 
most appropriate for ESR spectra of organic TFTs for the 



following reasons. First, the channel materials are com- 
posed of regularly aligned organic molecules that involve 
multiple degrees of freedom for nuclear spin moments. 
This feature clearly justifies our basic assumption that 
a single type of trap gives the Gaussian lineshape of the 
ESR spectrum. Secondly, it is possible to measure the 
high-precision ESR spectrum because of the fairly small 
spin-orbit interactions of organic materials. The field- 
effect device structure also enables the control of the car- 
rier density without introducing any randomness in the 
channel semiconductors. 

Such a direct probe is quite unique in investigating 
the microscopic carrier dynamics in the organic TFTs 
that have attracted considerable recent attention in the 
field of organic electronics. We have shown that the trap 
states in pentacene TFTs can be classified into three ma- 
jor groups: deep trap states with a spatial extension of 
about 1.5 molecules (A), relatively shallow trap states 
that extend over about 5 molecules (B), and shallower 
trap states that extend over 6 to 20 molecules (C). These 
states respectively correspond to deep and shallow trap 
states with binding energies of 140 meV (A) and 22 meV 
(B), with the broad feature ranging from 5 to 15 meV (C). 
These shallow in-gap states are crucial for understanding 
and improving the device performance of organic TFTs. 
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